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ABSTRACT 

We study the frequency distributions of S Scuti stars observed by the Kepler satel¬ 
lite in short-cadence mode. To minimize errors in the estimated stellar parameters, 
we divided the instability strip into ten regions and determined the mean frequency 
distribution in each region. We confirm that the presence of low frequencies is a prop¬ 
erty of all S Scuti stars, rendering meaningless the concept of S Sct/y Dor hybrids. 
We obtained the true distribution of equatorial rotational velocities in each region 
and calculated the frequency distributions predicted by pulsation models, taking into 
account rotational splitting of the frequencies. We confirm that rotation cannot ac¬ 
count for the presence of low frequencies. We calculated a large variety of standard 
pulsation models with different metal and helium abundances, but were unable to ob¬ 
tain unstable low-frequency modes driven by the k mechanism in any model. We also 
constructed models with modified opacities in the envelope. Increasing the opacity at 
a temperature logT = 5.06 by a factor of two does lead to instability of low-degree 
modes at low frequencies, but also decreases the frequency range of S Set-type pulsa¬ 
tions to some extent. We also re-affirm the fact that less than half of the stars in the 
S Set instability strip have pulsations detectable by Kepler. We also point out the huge 
variety of frequency patterns in stars with roughly similar parameters, suggesting that 
nonlinearity is an important factor in S Set pulsations. 
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1 INTRODUCTION 

Delta Scuti stars are main sequence dwarfs and giants with 
spectral types A0-F5 which pulsate in multiple p, g and 
mixed modes. They lie on the extension of the Cepheid in¬ 
stability strip towards low luminosities. From ground-based 
observations, most S Set stars are found to pulsate with fre¬ 
quencies in the range 5-50 d~^. The pulsations are driven 
by the k mechanism due to partial ionization of He II. 

The 7 Doradus stars lie in a fairly small region on, or 
just above, the main sequence that partly overlaps the cool 
edge of the <5 Set instability strip. They pulsate in multiple 
g modes with low frequencies (typically less than 5 d“^) and 
are thought to be driven by a mechanism known as “c on- 
vective blocking” JCuzik et al.ll2Oo'0 : lDupret et al.ll2004l l. A 
subset of the J Set stars pulsate in both the high and low- 
frequency ranges and are known as <5 Set/y Dor hybrids. 

Space photometry from the Kepler observatory has 
completely changed this view. It is found that all S Set stars 
pulsate in low -frequency modes, i.e. they are all hybrids, 
I^Balonall2014a^ . The misconception of hybrid <5 Set/y Dor 
stars is due to the fact that low frequencies in S Set stars have 


sufficiently high amplitudes to be observed from the ground 
only in a narrow band close to the red edge of the instabil¬ 
ity strip. The low frequencies in <5 Set stars present a serious 
problem because standard models predict that low frequen¬ 
cies are stable for S Set stars hotter than the granulation 
boundary. In these hot stars the subsurface convection zone 
is too thin for convective blocking to drive low-frequency 
pulsations. 


Apart from the problem of low frequencies in hot 
S Set stars, there is the additional problem that the ma- 
j ority of stars in the S Set in stability strip do not pulsate 
llBalona fc Dziembowskill2011^ . It is difficult to understand 
why in two stars with the same physical parameters one 
star should pulsate as a J Set star and the other not pul¬ 
sate at all, or at least have pulsational amplitudes not de¬ 
tectable in Kepler observations. It is possible that nonlinear 
mode coupling may stabilize the pulsations liDziembowskil 


_ i puisatior 

1982 :lDzie~mbowski~fc Krolikowskalfl9^:lDziembowski et al.l 


1988r) or that the opacity driving mechanism may be satu¬ 
rated dNowakowskil 2005l l . Since nonlinear models of nona- 


diabatic nonradial pulsations do not exist, it is not possible, 
at present, to test these ideas. 
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Another issue is the fact that A stars appear to be 
not as simple as previously thought. It has always been 
supposed that stellar activity ceases for stars hotter than 
the granulation boundary because the subsurface convec¬ 
tion zone is too thin to generate a magnetic field by the 
dynamo mechanism. Analysis of Kepler photometry shows 
that about 40 per cent of A stars have light variations whose 
periods closely match the expected rotation p eriods of these 
stars , suggesting the presence of starspots (jBalonal I2OIII . 
l2013h. Fnrthermore , about 2 per cent of A stars have flares 
( Balonall2012l . [201^ . In fact, the relative number of A-type 
flare stars is about the same as F and G flare stars an d not 
much smaller than K and M flare stars llBalonall2015l l. The 
starspot and flare activity in A stars indicates that we do not 
fully understand the physics of stellar envelopes which may 
be related to our lack of understanding of low frequencies in 
hot S Set stars. _ 

Using Kepler data. iBalonal (l2014al l showed that low fre¬ 
quencies in S Set stars are not stochastically driven, nor 
can they be explained as nonlinear combinations of high- 
frequency modes. By comparing the observed distribution 
of frequencies i n Kepler S Set stars with those from models, 
IBalonal (l2014al l also showed that the low frequencies can¬ 
not be explained as a result of rotational splitting of high- 
frequency modes. 

In order to understand the origin of the low frequen¬ 
cies in <5 Set stars, it is necessary to explore what effects 
contribute most to the stability or instability of these fre¬ 
quencies in the models. A simple examination of the insta¬ 
bility parameter, rj, calculated by a nonadiabatic pulsation 
model will tell us if a particular mode is stable or unsta¬ 
ble. However, in order to compare the predicted frequency 
distribution with observations, it is necessary to simulate ro¬ 
tational splitting of the frequencies as well as other effects. 
In the absence of mode identification and sufficiently accu¬ 
rate stellar parameters, this is the only way to test if model 
predictions agree with observations. In order to simulate the 
effect of rotational splitting, the distribution of true equato¬ 
rial rotational velocities for t he group of stars under study 
is required. In lBalon'^ ll2014ah all S Set stars were treated as 
a single group using rotation periods inferred from the light 
curve. 

In this paper we compare the observed and simulated 
frequency distributions in several regions across the insta¬ 
bility strip. Different physical conditions are used for each 
set of models in an attempt to determine which effects most 
closely reproduce the observed frequency distributions. In 
this way we hope to determine whether it is possible to un¬ 
derstand the low frequencies under the assumption that they 
are driven by the k mechanism. We study mode stability 
using standard models, i.e. using different abundances and 
standard opacities, as well as non-standard models, i.e. using 
opacities which are artihcially increased in certain regions in 
the envelope. 


2 THE DATA 

The Kepler light curves are available as uncorrected simple 
aperture photometry (SAP) and with pre-search data con¬ 
ditioning (PDC) in which instrumental effects are removed. 
The vast majority of the stars are observed in long-cadence 


(LC) mode with exposure times of about 30 min. Several 
thousand stars were also observed in short-cadence (SC) 
mode with exposure times of about 1 min. These data are 
publicly available on the Barbara A. Mikulski Archive for 
Space Telescopes (MAST, archive.stsci.edu). 

In order to identify <5 Set stars in the Kepler held, we 
calculated the periodograms of over 20 000 stars, including 
all stars observed in SC mode. By visually examining the 
light curves and periodograms and using the effective tem¬ 
perat ures in the Kepler Input Catalogue IKIC. lBrown et ^ 
I 2 OIIII as a guide, we were able to identify over 1600 S Set 
stars in LC mode and 403 S Set stars in SC mode with 
known stellar parameters. These do not include 5 Set stars 
in eclipsing binary systems. 

The highest pulsation frequency that can be detected in 
LC mode is about 24 d~^. Since the pulsation frequencies in 
many 5 Set stars exceed this value, only SC observations al¬ 
low the unambiguous identiheation of all pulsation frequen¬ 
cies. Due to variations in heliocentric time correction, the 
Kepler data are not sampled at exactly equal time intervals. 
In principle, this allows frequencies higher th an 24 d“^ to 
be identihed in LC data (iMurphv et al. l2013ll . In practice, 
however, the chance of frequency misidentiheation is very 
high for low-amplitude peaks because the difference in peak 
amplitude between the high- and low-freqnency aliases is 
slight. For this reason, we restrict our analysis to the 403 
S Set stars observed in SC mode. 

To extract all significant frequencies we hrst calcu¬ 
lated the standard Lomb per iodogram for uneqnally-spaced 
data dPress &i Rvbickj Il989ll . The usual technique to ex¬ 
tract peaks is that of successive prewhitening. Great caution 
needs to be exercised in using this method as indiscr iminate 
use lea ds to a large number of spurious frequencies (IBalonal 
l2014bll . In effect, only those peaks which are visible in the 
periodogram of the raw, un-prewhitened data can be consid¬ 
ered as significant. Indiscriminate successive prewhitening 
will extract frequencies which cannot be resolved in the peri¬ 
odogram, even though the extracted amplitudes of such spu¬ 
rious frequencies may be deemed significant. In determining 
the frequency content of a star, we terminated prewhiten¬ 
ing when the peak can no longer be resolved in the original 
periodogram. 


3 SYSTEMATIC AND RANDOM ERRORS 


In order to compare the observed frequency distribution 
with the calculated frequency distribution we need to de¬ 
termine the accuracy of the stellar parameters for individual 
stars. The stellar parameters derived from Sloan multicolour 
photometry (without the u band) are listed in the KIC. It is 
important to determine the systematic and random errors 
that might be present in the values of effective tempera¬ 
ture, Tefj, and relative luminosity, logL/LQ, derived from 
the KIC parameters. 

We do not know the true values of Tefi and logL/L© 
for any 5 Set star, but we can compare the values in 
the KIC derived from multicolour photometry with the 
more prec ise values that can be d erived from spect r oscopy . 


more prec ise values tnat can be d erived trom spect r oscopy . 
Recently, Lehmann et ahl (I2OIIII: Catanzaro et al.l (1201(31 : 
iNiemezura et all ( 2015l l and Tto,chenko et al.l ( 2013fl have 
obtained spectroscopy of several Kepler stars from which 
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Figure 1. Top panel: the corrected KIC effective temperature 
as a function of the effec tive temperature determ ined by spec¬ 
troscopic observations of iNiemczura et aP ll201fi1 . The straight 
line has unit slope and zero intercept. Bottom panel: the value 
of log(L/L 0 ) derived from the KIC parameters as a function 
of log(L/LQ) derived from the spectroscopic parameters. The 
straight line has unit slope and zero intercept. 



Figure 2. The theoretical H-R diagram showing the SC 5 Set 
stars (dots) and the designated regions (labeled). The line is the 
theoretical zero-age main sequence. 


able, the luminosity, log(L/Z/ 0 )(KIC) can likewise be esti¬ 
mated from the tabulated effective temperature and radius. 

The bottom panel of Fig. [T] shows log(L/L 0 )(spec) as 
a function of log(Z//Z/ 0 )(KIC). The correlation is rather 
poor, no doubt due to the poor precision of lo g(ff) from 
the photometry. According to iBrown et al.l (1201 ll) the stan¬ 
dard deviation of log g in the KIC for dwarfs is about 
0.4 dex, whereas it is a bout 0.15 dex from the spectroscopy 
llNiemczura et al ][ 2 oTi). As the figure shows, a line of unit 
slope and zero intercept gives a satisfactory fit. The stan¬ 
dard deviation in \og{L/L q) is 0.26 if the regression is per¬ 
formed with log(L/L 0 )(KIC) as the independent variable or 
0.29 as the dependent variable. The standard deviation of 
\og{L/ Lq) derived from the KIC parameters must therefore 
be about 0.20 dex, since the total variance is the quadratic 
sum of the two variances. 

We used spectroscopically derived parameters whenever 
possible for each star. When only KIC values are known, we 
corrected Tefj listed by the KIC by 144 K and estimated the 
luminosity from the KIC radius and uncorrected effective 
temperature. 


stellar parameters a re derived. The most nu merous obser¬ 
vations are those of INiemczura et all (( 201 ^ for which we 
find a systematic difference in effective temperature ATefr = 
Teff(spec) — reff(KIC) = 144 ± 24K from 107 stars with a 
standard deviation of 253 K per star. There is no significant 
dependence of ATefi o n Teg. Since the stand ard deviation 
of Teff(spec) « 150 K (INiemczura et al.|l2015l ) , we estimate 
that the standard deviation of Teg (KIC) « 200 K on the 
assumption that the variances add in quadrature. The top 
panel of Fig.[T]shows reff(spec) from lNiemczura et al.l (l2015l ) 
as a function of reff(cor) = Tefr(KIC) -I- 144. 

Given Teg and the surface gravity, \ogg, t he stellar ra¬ 
dius m ay be estimated using the relationship in I Torres et al.l 
(l2010t) . This relationship requires the metal abundance 
which is usually listed among the spectroscopically derived 
parameters. When the metal abundance is not available, the 
solar value was used. From the effective temperature and the 
radius, the relative stellar luminosity, log(L/I/ 0 )(spec), can 
be determined. For stars where only KIC values are avail¬ 


Given the large errors in Teff and \ogL/LQ, the lack 
of mode identiheation, and the profound effect of rotation 
on the frequency spectrum, it is not possible to model the 
pulsation frequencies of individual stars. In fact, we are not 
interested in detailed modeling of a specific star, but only 
in a comparison of the overall observed and calculated fre¬ 
quency distributions. For this reason we have taken a dif¬ 
ferent approach. We can minimize the uncertainties in effec¬ 
tive temperature and luminosity by using the mean values 
of log Teff and \ogL/LQ of many stars in approximately the 
same location in the instability strip. A good estimate of the 
mean frequency distribution at this position in the instabil¬ 
ity strip is obtained by averaging the frequency distributions 
of these stars. For this purpose we divided the theoretical 
H-R diagram into ten regions as shown in Fig. [21 These re¬ 
gions were chosen to cover the instability strip as evenly as 
possible with reasonable resolution while, at the same time, 
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Table 1. For each region number, the number of S Set stars, N, in that region and the average number of modes, Af^iodesi 

having u > 5d“^ are shown. Also shown is the corresponding range in spectral type, the mean effective temperature, Tefi, luminosity, 
log L/Lq, surface gravity, logg and radius, R/Rq. The last column, (ue) is the mean equatorial rotational velocity (in kms“^) derived 
from observations of many field stars of the particular spectral type and luminosity class range. 


N-Reg 

N 

-^modes 

Sp.Ty. 

log Teff 

logL/Le 

logs 

R/Rq 

{Ve) 

1 

29 

90 

F8V-F3V 

3.8253 + 0.0015 

0.6889 + 0.0152 

4.100 + 0.011 

1.654 + 0.024 

40 

2 

60 

252 

F4V-A9V 

3.8546 ± 0.0012 

0.8896 + 0.0119 

4.043 + 0.009 

1.822 + 0.020 

90 

3 

32 

183 

F0V-A7V 

3.8858 + 0.0019 

1.0922 + 0.0178 

4.003 + 0.011 

1.993 + 0.030 

130 

4 

25 

138 

A8V-A4V 

3.9146 + 0.0012 

1.2671 + 0.0190 

3.985 + 0.013 

2.136 + 0.039 

150 

5 

35 

238 

F5IV-F1IV 

3.8537 + 0.0014 

1.0808 + 0.0136 

3.869 + 0.010 

2.280 + 0.030 

70 

6 

61 

245 

F1IV-A8IV 

3.8837 + 0.0012 

1.3338 + 0.0130 

3.781 + 0.010 

2.663 + 0.035 

110 

7 

34 

220 

A8IV-A5IV 

3.9121 + 0.0014 

1.4897 + 0.0148 

3.777 + 0.013 

2.793 + 0.047 

140 

8 

18 

103 

F5III-F1III 

3.8471 + 0.0022 

1.3958 + 0.0210 

3.565 + 0.016 

3.381 + 0.070 

80 

9 

58 

365 

F1III-A8III 

3.8746 + 0.0013 

1.5508 + 0.0118 

3.551 + 0.009 

3.560 + 0.041 

100 

10 

15 

83 

A9III-A6III 

3.9043 + 0.0025 

1.7030 + 0.0181 

3.553 + 0.012 

3.692 + 0.059 

110 


allowing a fairly large number of stars in each region. The 
mean effective temperature, luminosity, surface gravity and 
radius of all <5 Set stars in each region are shown in Table[T] 


4 MODE VISIBILITY 

In determining the frequency distribution of a particular star 
we are faced with the problem that individual frequencies are 
of widely different amplitudes. It would, of course, be prefer¬ 
able to calculate the amplitudes from the models as well. 
Unfortunately nonlinear calculations of nonradial modes are 
beyond existing technical capabilities. One approach is to 
take the linear growth rate as a proxy for the amplitude. 
However, the final pulsational amplitude attained by a star 
depends on many factors and not only the initial growth 
rate. We decided that the most satisfactory method is to 
ignore amplitude altogether. In other words, a mode of very 
high amplitude is given the same weight as a mode of very 
low amplitude. In determining the observed frequency distri¬ 
bution, we simply count the number of significant frequency 
peaks in the periodogram within a given frequency interval, 
ignoring their amplitudes. 

Using this approach requires that each unstable mode 
in the models be weighted according to its visibility. It is 
evident that modes of very high spherical harmonic degree, 
I, will not be observed because of cancellation effects and 
cannot therefore be given the same weight as modes of low 
1. In other words, each unstable mode derived from a model 
should be assigned a weight equal to its visibility before it 
can be compared with observations. To calculate the visi¬ 
bility of a mode, we need to discuss the expected relative 
light amplitude caused by a mode of particular spherical 
harmonic degree, I, and azimuthal order, m. 

The monochromatic light amplitude, A\{i), for a pul¬ 
sating star whose axis of rotation is inclined at an- 
gle i with respect to the observe r can be expressed as 
IIDaszvhska-Daszkiewicz et al.ll2002l l: 

Axii) = eVrii, 0)btRe{Dli + Da,, -f D^,,}. 

Here e is the intrinsic amplitude, Y{^(i,(f)) = 

IVimD,’”(cos is the spherical harmonic with 

the normalizing factor and P,™'(cosi) the associated Legen¬ 
dre polynomial (i = x/—l). The disc averaging factor, 


can be calculated for any particular limb-darkening law. In 
the above expression, the Di,, term describes the effect of 
temperature variation on the light amplitude, whereas the 
influence of effective gravity changes is contained in the 
D 3 ,, term. The D 2 ,, term describes the effect of variations 
in projected area on the light amplitude. These terms are 
given by 


X 1 d\ozTx\hn 

4-^ aiog(reff) 


D2,i = ( 2 + 0 ( 1 - 0 , 


D 


A 

3, 


I — — 



dlogTx\bt\ 

i91og(c,efr) 


The ratio of luminous flux to displacement, /, for each mode 
is calculated by the nonadiabatic pulsation code. The partial 
derivatives of the luminous flux, Rx, with respect to effec¬ 
tive temperature, Teft, and with respect to effective gravity, 
geH, can be obtained from model atmospheres. The angular 
pulsation frequency is uj, and p is the mean density of the 
star. 

If we consider a group of stars in the same region, the 
visibility of a particular mode for the group is obtained by 
averaging the various terms in the above equation. Because 
we are only counting unstable modes, irrespective of ampli¬ 
tude, e is unity if the mode is unstable and zero if it is stable. 
The average value of 0) is found assuming random ori¬ 
entation of the axis of rotation so that, putting x = cos(z). 


with 


(yr(*,o)) 


= ivr(|p,"‘(cosi)i), 

_ 2 A \prix)\ 


dx, 


Nr = 


{2l + l){l-my. 
4 :Tt{ 1 + m)\ 


The value of (Y,"‘(i,0)) does not change very much and lies 
within the range 0.13-0.30 for modes with I + 6 . 

The disk-averaging factor, , is given by 

bi= hx{n)iJ-Pi{n)dn, 

Jo 

where g, = cos 6 and 6 is the angle in the spherical coordi¬ 
nate system centered on the star. We used the simple limb- 
darkening law h{g) = 1 + |/i independent of wavelength. 
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Values of bi drop sharply with Z: bo = 1, &i = 0.708 ,.. .bo = 
0.008. 

In calculating Di^i and we used the par- 

tial derivatives for the Johnson V band computed by 
iKowalczuk &i Daszvhska-Daszkiewicd (l2007h . The value of 
Di,i is typically around 40 and is in the range 2-3. The 
geometric term, D 2 ,i increases in absolute value with I and is 
comparable with Di^i for high 1. The temperature variation 
together with bi dominates the visibility. 


5 THE EFFECT OF ROTATION 

In a non-rotating star, pulsation modes with the same spher¬ 
ical harmonic degree, Z, but with different azimuthal num¬ 
bers, m, have the same frequency. To first order, the effect 
of rotation is to remove this degeneracy so that a mode 
with frequency vq in the non-rotating star appears as 2Z -|-1 
equally-spaced multiplets with frequencies, v, given by 

U = 120+ mUrot (1 - Cnl) , 

where Vrot is the rotation frequency and Cni is a constant 
which depen ds on the stru cture of the star, Z and the ra¬ 
dial order, n (lLedouxlll95ll ~). The value of Cni is quite small 
for p modes and hence the frequency splitting may be sev¬ 
eral cycles per day in the extreme case of high m and large 
r'rot. The effect of rotation to first order is to introduce a 
symmetric spread in frequencies. 

The above formula breaks down as the rotation rate 
increases and the splitting is no longer symmetric. The fre¬ 
quency of each multiplet decreases with increasing rotation 
rate, but the decrease in frequency is greater for larger values 
of |m|. Thus the spread in frequencies is skewed towards low 
frequencies. It is therefore important to take this effect into 
account by including higher- order rotation al perturbations. 

The following formula, jGoupill I 2 OIIII , is accurate to 
third order: 

to = LUO — m(l — Cni)+l + {Di + rrd + m{Ti + npT 2 )fl^, 

where to, uo are the perturbed and unperturbed angular pul¬ 
sation frequencies and Q is the angular rotation frequency. 
The values of Di, D 2 , Ti, and T 2 depend on the mode. In 
this expression the pulsation and rotation frequencies are in 
units of the dynamical frequency Wdyn = l/GM/B+ where G 
is the gravitational constant, M the stellar mass and R the 
stellar polar radius. The coefficients have not bee n calculated 
for sta rs in the <5 Set instability strip. However. iReese et al.l 
ll2006ll lists values for a polytrope of index 3 for Z Si 3 and 
n ^ 10. 

Since most A stars are moderate or rapid rotators, it 
is very important to include rotation effects as accurately 
as possible. Since calculations are only available for poly¬ 
tropes and for limited values of Z and n, we are faced with 
a problem. We decided that the best approach is to use the 
third-order formula and polytropic values for all p modes. 
For p modes with Z > 3 we use the coefficients for Z = 3. For 
radial orders greater n > 10 we use the values for n = 10 . 
This might overestimate the rotational splitting for modes 
with Z > 3. Rotational splitting for g modes is smaller than 
for p modes. We decided to use the first-order formula for 
all g modes. The distinction between p and g modes is made 
on the basis of the ratio of the mode kinetic energy in the 


gravity-wave propagation zone, Etg, to the total kinetic en¬ 
ergy of a mode, Ek, as calculated by the models. We treated 
modes with Ekg/Ek < 0.5 as p modes. 

To calculate the frequencies of the rotational multiplets 
in the way just descri bed requires kn owledge of the rota¬ 
tion frequency, r'rot. In lBaloria (l2014al ) the value of r'rot was 
obtained from a peak in the periodogram which could be 
attributed to rotational modulation. We cannot do this be¬ 
cause the number of stars with known photometric rotation 
period is too small (and sometimes zero) in each region. 
To overcome this problem we have to make the assumption 
that the distribution of equatorial rotational velocities in a 
particular region closely corresponds to the distribution of 
equatorial rotational velocities for field stars with a spectral 
type and luminosity class appropriate to that region. We 
used the catalogue of projected rotatio nal velocities com¬ 
piled by iGlebocki fc Stawikowskil ll2000ll to determine the 
distribution of v sin i for stars in each region. To obtain the 
distribution of true equato rial rotational velocities, we used 
the procedure described in lBalon'3 (Il975l l. Fig.[3] shows the 
distribution of usini and the polynomial approximation to 
the true distribution of equatorial rotational velocities. 


6 ASSIGNING WEIGHTS 

For a stellar model with a given mass and radius we may 
calculate the rotational frequency, iZrot, for any given equa¬ 
torial rotational velocity, v^. Given i/rot and the calculated 
frequency of a unstable mode of degree Z and radial order 
n, we may calculate the 2Z + 1 rotationally-perturbed fre¬ 
quencies using the first-order Ledoux formula for g modes 
and the third-order formula for the p modes. To each of 
these frequency multiplets we assign the same weight, Wrot, 
which is the probability of finding a star with the particular 
value of Ve- This weight is calculated using the polynomial 
fit to the distribution of true equatorial velocities for stars 
in the appropriate region. We calculated the frequencies and 
weights of all rotationally-split multiplets using values of 
from zero to the maximum value in steps of lOkms”^. 

We also assign a weight Wi = (Yi’^(i, 0)) which accounts 
for the visibility of the mode due to random orientation of 
the axis of rotation by numerically calculating the integral. 
This term is relatively unimportant because the Wi does 
not change much. The visibility due to cancellation effects, 
Wb = bi, is pre-calculated for any given value of Z. Finally 
we calculate the weight Wf = + D 2 ,i + Do,i using the 

values of f,LU and p from the model as well as the partial 
derivatives applicable to the given model. The total weight 
for the mode with frequency 12 is Wu = WrotWiWbWf. In 
calculating the predicted frequency distribution we chose a 
bin size of ld“^. The probability at any given frequency, 
P{i 2 ), is P{i 2 ) — X) Wi, for all frequencies within this interval. 

In principle, the frequency distribution, Piy), needs to 
be calculated for all possible values of 1. However, the mode 
visibility drops sharply with Z, so that modes with Z > Zmax 
have amplitudes which are so low that they can be ignored. 
The simplest way to determine Zmax is to construct frequency 
distributions using increasing values of Zmax- We find that 
there is scarcely any difference in the distributions with 
Zmax > 4. We adopted Zmax = 6 which is certainly adequate 
for our purposes. 
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Figure 3. The distribution of projected rotation velocity, vsini, 
for field stars with spectral types appropriate for each region is 
shown by the histogram. The distribution of true equatorial ro¬ 
tational velocities, Ue, is shown by the dashed curve. 


In the above procedure we have made assumptions that 
may not perhaps be entirely justifiable. For example, we 
assume that all rotationally-split multiplets have the same 
intrinsic amplitude. We cannot calculate these amplitudes 
so we do not really know, but this is probably a fair as¬ 
sumption when averaged over many stars. It seem reasonable 
to assume that stars with approximately the same stellar 
parameters and rotational velocities have similar frequency 
distributions. This is certainly the case given our current un¬ 
derstanding and that is what the models predict. We shall 
see below that even this seemingly secure assumption may 
not be correct. 


7 COMPARISON WITH THEORETICAL 

PREDICTIONS 

Equilibrium models w ere computed us ing the Warsaw - New 
Jersey evolution code (|Paczvnskilll97d l. assuming two initial 
hydrogen fractions, Xo = 0.70, 0.65, and four metal abun¬ 
dances, Z = 0.005, 0.015, 0.030 and 0.050. We adopted the 
chemical element mixture of lAsplu nd et al] ||2009|) and two 
sources of the opacity data: OPAL ( Rogers fc Iglesi^ll992ll 
and OP llSeatorj2005lL We did not include overshooting from 
the convective core and used the mixing length parameter 
Q!mlt = 0.5 for the convective scale height. We chose this 
low value of q:mlt to eliminate, as far as possible, the fic¬ 
titious instability of low-frequency modes in the cooler and 
lower-mass models. This instability is a consequence of the 
frozen convective flux approximation that we adopted. 

We calculated a grid of equilibrium models covering the 
S Set instability strip. The models take rotation into ac¬ 
count to first order by applying the centrifugal force cor¬ 
rection to local gravity while keeping spherical symmetry. 
We used a typical initial equatorial rotational velocity of 
Ve = 150kms“^ except for the models with the lowest 
mass where we used Ue = 50kms~^. The effect of rota¬ 
tion is a slight shift of the evolutionary tracks to lower ef¬ 
fective temperatures and higher luminosities. The main se- 
quence band is also somewh at extended (see, for example, 
iBreger fc Pamvatnvkhiri998l L Given the fact that we need 
to determine the rotational frequency splitting over a wide 
range of rotational velocities corresponding to the known 
equatorial rotational velocity distribution for a particular 
group of stars, the value of Vs in the equilibrium model is 
not an important factor. A non-rotating model would serve 
our purposes equally well. As explained in the previous sec¬ 
tion, we used the 3-rd order rotational perturbation only for 
p modes, while for g modes the 1-st order Ledoux splitting 
formula was used. _ _ 

The linear nonadiabatic code of iDziembowskil (Il977l j 
was used to determine the frequencies and the instability 
parameter, rj, as well as the the ratio of luminous flux to 
displacement, /, for each mode. Only modes with I 6 
were considered. 

In calculating the predicted frequency distribution, we 
use individual values of log Teg and log L/Lq for each star 
in a particular region. The model with temperature and lu¬ 
minosity closest to these values is used to determine the pre¬ 
dicted frequency distribution. The average of all frequency 
distributions for stars in a given region is taken as the best 
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Figure 4. The observed frequency distributions (filled) and the corresponding theoretical frequency distributions derived from pulsation 
models with metal abundance Z = 0.015 (solid blue curves) and Z = 0.030 (dotted green curves). Only unstable modes with the degrees, 
I ^ 6 are considered. Each panel shows the frequency distribution corresponding to the labeled region in the H-R diagram of Fig. [2] 


approximation of the predicted frequency distribution for 
that region. 

The observed and predicted frequency distributions are 
shown in Fig.|4] for Z = 0.015 and Z = 0.030 using OPAL 
opacities and an initial hydrogen abundance Xq = 0.70. The 
unstable low frequencies predicted for the models of region 1 
are a result of the frozen convection approximation and are 
fictitious. It can be seen that low frequencies are observed 
in stars across the whole instability strip. In other words, 
all (5 Set stars are hybrids. Also, it can be seen that the 
predicted frequency distributions fail to show frequencies 
less than about 5d“^. The same result was obtained with 
models calculated using OP opacities. 

Apart from the disagreement in the low-frequency 
range, we also note that for the more luminous stars of re¬ 
gions 8, 9 and 10 the observed distributions extend to much 
higher frequencies than predicted. In general, we expect the 
observed distribution to be broader than the predicted dis¬ 
tribution because each region includes a fraction of stars 
with true values of Tefi and logL/I/Q which are outside the 


region. Even so, it does not seem possible to explain the 
high-frequency tails of the luminous stars in this way be¬ 
cause there is no such tail among the less luminous stars. 

What we can definitely conclude from Fig.H] is that a 
change in metal abundance cannot explain the mismatch 
between observations and the models. We also calculated 
models with Z = 0.005 and Z = 0.050 which we do not 
show in Fig.[4] to avoid confusion. The distributions with 
lower and higher metal abundance are quite similar to those 
shown in the figure. We also calculated models where the 
helium abundance is substantially increased so that the ini¬ 
tial hydrogen abundance is decreased from Xq = 0.70 to 
Xq = 0.65, but this does not affect mode stability at all at 
low frequencies. However, there is a slight increase of the in¬ 
stability parameter, 77 , at high frequencies. We also studied 
the effect of modes of higher degree and found a negligible 
difference in the distributions for modes with Zmax = 6 and 

^max — 10. 

It could be argued that we have not taken the effect 
of rotation fully into account and that agreement could be 
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Figure 5. The same as Fig.|4] but we have artificially increased the effect of rotational splitting by using the 3-rd order formula for 
all modes and setting the disc-averaging factor, fe;, to unity so as to maximize the visibility of modes of high degree where rotational 
splitting is greatest. Only predicted distributions for Z = 0.015 are shown. 


achieved by a more realistic treatment of rotational split¬ 
ting. Unfortunately, we do not have access to non-adiabatic 
models which take into account stellar distortion, gravity 
darkening and rapid rotation in a more realistic way. We 
can, however, take the effect of rotation to extremes by us¬ 
ing the 3-rd order perturbation for all modes. We can also 
enhance the visibility of modes of high degree (which have 
maximum rotational splitting) by artificially suppressing the 
disc-averaging factor, bi (i.e, setting 6; = 1 for all modes). 
The distributions calculated in this way are shown in Fig. [5] 
for Z = 0.015. 

Although there is a slight increase in the numbers of 
low-frequency modes, it is still far too little to match obser¬ 
vations. Artificially increasing the effect of rotation in this 
way also does not assist in explaining the high-frequency tail 
in luminous stars. We may confidently conclude that there 
is no possibility that rotational splitting of high-frequency 
modes can can explain the low frequencies. 

Since the low frequencies cannot be explained by modi¬ 


fications of standard models, we need to consider other pos¬ 
sibilities. 

Low frequencies not predicted by models also occur in 
the /? Cep stars. These are early B stars with multiple p, 
g and mixed modes driven by the opacity bump due to 
iron group elements (the Z bump). Some low frequencies 
in these stars can be explained by standard models and 
some demand increasing the opacity in the Z-bump region 
llPamvatnvkh et al.ll2004h . Whether this is the correct solu¬ 
tion or how to create the increased opacity is not known. As 
in the B stars, it is possible that artificially increasing the 
opacity of S Set models in the Z-bump region at log T « 5.35 
might destabilize the low frequencies. 

Recently, ICueie 3 (|2012| . |2014| ~) found that the OPAL 
and OP opacities are markedly underestimated in com- 
parison with the Rossela nd mean opacities taken from the 
ICastelli fc Kuru^ (l2003li model atmospheres. As a result, a 
new opacity bump appears at logT « 5.06. This bump is 
due to an att empt to inclu de all spectral lines of atoms in 
the opacities (lKurucjl201lll . This additional opacity bump 
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Figure 6. The instability parameter, rj, as a function of frequency 
for modes with Z ^ 6 for representative models of region 3. All 
models have logTefj = 3.8855 and logL/LQ = 1.093. Panel (a) is 
a standard model with mass M = 1.76 Mq and Z = 0.015. Panel 

(b) is a standard model with M = 2.25 Mq and Z = 0.050. Panel 

(c) is a standard model with M = 1.48 Mq and Z = 0.005. In 
panel (d) the opacity at logT = 5.35 (the Z-bump region) has 
been increased by a factor of two. In panel (e) the opacity at 
logP = 5.06 has been increased by a factor of two to simulate a 
new opacity bump in the Kurucz opacity data (see main text). 
All models have initial hydrogen abundance Xq = 0.70. 


affects the stability of stellar pulsations and also needs to be 
considered. In testing the effect of this bump (which we will 
call the “Kurucz” bump), we did not use the actual Kurucz 
opacities where the bump occurs, but simulate this addi¬ 
tional opacity bump by artificially increasing the standard 
OPAL opacities in the appropriate temperature range. 

In Fig.[S]we show the instability parameter, rj, as a func¬ 
tion of frequency for Z ^ 6 in five models with the same effec¬ 
tive temperature and luminosity but different metal abun¬ 
dances and enhanced opacities. As one can see, an increase 
in 77 in the low frequency range can be obtained either by 
increasing the value of Z or by an increased opacity in the 
Z-bump or by increasing the standard OPAL opacities to 
simulate the Kurucz bump. However, low-frequency modes 
remain stable except when the Kurucz bump opacity is in¬ 
cluded. Instability of Z = 1 modes at low frequencies begins 
when the standard OPAL opacity in the Kurucz-bump re¬ 
gion is increased by a factor of two. Maximum instability 
occurs when this opacity is increased by a factor of three. 
Any further opacity increase results in saturation and no fur¬ 
ther increase in rj can occur. Note that an increase in opacity 
in the bump or Kurucz-bump regions leads to somewhat 
reduced rj for modes with high frequencies. For example, 
the frequency range of unstable high-frequency modes is de¬ 
creased from the 15-35 d“^ range for a standard model to 
16-32 d“^ for a model where the OPAL opacity in the Ku¬ 
rucz bump region is increased by a factor of two. 


8 FURTHER UNSOLVED PROBLEMS 

Even a cursory inspection of the periodograms of <5 Set stars 
is sufficient to show the wide variety of frequency patterns 
in these stars. In fact, each star is unique and can be identi¬ 
fied by its periodogram. This is strange because one expects 
stars with similar effective temperatures, luminosities and 
rotational velocities to have similar frequency spectra, but 
this does not seem to be the case. Examples of this effect 
are shown in Fig.[7]for stars belonging to region 3. The stars 
shown all have well-determined spectroscopic effective tem¬ 
peratures and surface gravities. 

While it is true that the stars may individually have 
considerably different parameters due to substantial errors 
in Teff and logL/Lg, the differences in frequency spectra 
are very pronounced. This disparity cannot be proved until 
more accurate parameters are obtained, of course. Neverthe¬ 
less, it seems that small differences in the parameters may 
lead to large differences in frequency patterns, implying that 
nonlinearities in S Set envelopes may be very important. 

The question of the number of non-pulsating 
stars in the S Set in s tabilit y strip was addressed by 
iBalona fc Dziembowskil (l201ll ). They found that most 
5 Set stars have effective temperatures in the range 
7000 < Teff < 8500 K and that even in this range no more 
than 40-50 per cent of stars pulsate as S Set variables. We 
can confirm this finding. We have found 1165 <5 Set stars 
in this temperature range in the Kepler field observed in 
long-cadence mode. There are 2839 stars (including the 
S Set stars) in the same temperature range, meaning that 
only about 41 per cent of stars in the instability strip are 
detected as <5 Set stars in the Kepler photometry. 

It could be argued that most of the non-5 Set stars 


© 2011 RAS, MNRAS OOO.fTHTTI 
































10 L.A. Balona and J. Daszynska-Daszkiewicz 


4 

2 

0 

3 

2 

1 

0 

2 

Q. 1 
Q. 

■Q 

i 0 

Q. 

< 

2 


0 
3 

2 

1 

0 

1 

0.5 
0 

0 5 10 15 20 25 30 35 

Frequency (d'^) 

Figure 7. Examples of periodograms of <5 Set stars with similar 
effective temperatures and luminosities (in region 3). All stars 
(KIC numbers are shown) have spectroscopic determinations of 
effective temperature and surface gravity. 



are outside the instability strip due to errors in the effec¬ 
tive temperature. Errors in luminosity alone play no role 
because the stars will still be within the instability strip as 
they can only lie between the ZAMS and TAMS. There are 
no known high-luminosity stars in this temperature range 
and in any case none are expected due to the very short 
lifetimes such stars will have in this evolutionary stage. We 
have seen that a typical error of about 200 K in Teff can 
be expected, which is considerably less than the range of 
1500 K that we are discussing. The probability that a star 
in the middle of the instability strip lies, in actual fact, out¬ 
side the strip is less than 0.01. Of course, the probability 
will be higher if the star is closer to the edge of the insta¬ 


bility strip, but it means that the probability that all 1674 
non-pulsating stars are outside the instability strip is the 
product of the individual probabilities which is essentially 
zero. There can be no question that non-pulsating stars are 
present in the 5 Set instability strip. Of course, it can be 
argued that the pulsations are less than the detectable limit 
for Kepler (typically less than about 50ppm). This might 
be the case, though one can argue against this idea on the 
basis of the amplitude distribution (iBalona &: Dziembowskil 
I 2 OIII) . Such an amplitude disparity, if it exists, still needs 
an explanation. 


9 CONCLUSIONS 


We have compared the frequency distributions of groups 
of S Set stars with similar effective temperatures and lu¬ 
minosities with the frequency distributions calculated from 
pulsation models taking into account, in an approximate 
but fairly realistic way, the effect of r otational splittin g. We 
have come to the same conclusion as iBalonal (l2014al ) : that 
it is impossible to account for frequencies with u < 5d~^ 
on the basis that these are high-frequency modes shifted to 
low fre quencie s by rot ation. Other factors examined in de¬ 
tail in IBalonal ll2014ah are also excluded. The presence of 
low frequencies is a general feature of all <5 Set stars, no 
matter where they lie in the instability strip. This is clearly 
seen in the observed distributions shown in Fig.|4l We also 
found other inconsistencies between the observed and pre¬ 
dicted frequency distributions. For example, the observed 
distribution for the more luminous stars is much wider than 
expected. However, this is a minor problem compared to 
absence of predicted low frequencies. 

We calculated a large variety of standard models with 
differing helium and metal abundances, but in no case could 
we find unstable modes at low frequencies. The instability 
parameter, rj, does however tend to increase as the metal 
abundance, Z, increases. We also studied the effect of in¬ 
creasing the opacity at temperatures logT = 5.35 (the Z 
bump) and the effect of simulating a new opacity bump at 
log T = 5.06 (the Kuruez bump). The Kuruez bump does not 
occur in the OPAL and OP opacities, but appea rs in model 
atmospheres, as discussed bv ICugieil 1 20121 . l2014h . We found 
that increasing the opacity in these two regions does increase 
the value of rj. However, low-frequency modes remain stable 
unless the OPAL opacities at logT = 5.06 are increased by 
at least a factor of two to simulate the Kuruez bump or if 
the opacity in the Z bump is increased by a factor of at least 
three. At the same time, the range of high-frequency modes 
is decreased to some extent. We do not claim that this opac¬ 
ity increase is the solution to the problem of low frequencies 
in S Scuti stars, but only that there is a likely problem with 
current opacity data and that further investigation into the 
sources of opacity in this region are required. 

A problem of equal importance is the question of why 
there are so many non-pulsating stars in the 5 Set instability 
strip. This could simply be that the Kepler photometry is 
not sufficiently precise to detect such pulsations. This does 
not resolve the problem because it would not explain the 
very large disparity in amplitudes among the 5 Set stars. We 
are clearly faced with substantial problems in the physics of 
A-F stars. 
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Recently, further problems in our current understanding 
of stellar pulsations have come to light. It seems that there 
are a group of pulsating stars with high frequencies lying 
between the cool end of the /3 Cep instability strip and the 
hot end of the S Set instabil it y stri p. These were recently 
investigated by iBalona et al.l ( 201!Th using Kepler and K2 
data. Standard models cannot reproduce these pulsations, 
though it is possible that these so-called “Maia” variables 
may be rapidly-rotating SPB stars. The unexplained pres¬ 
ence of low frequencies in a well-studied group such as the 
5 Set stars shows that we still do not fully understand the 
physics and/or envelopes of hot stars. Until we have a bet¬ 
ter understanding of the low frequencies in A stars we are 
not likely to make much progress with the more complex 
problem of explaining the origin of high frequencies in Maia 
variables. 
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